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^ ■ Abstract. Collective phenomena from RHIC and LHC experiments indicate a 

_ strongly coupled quark gluon plasma. Finite temperature lattice QCD calculations 

show the interaction measure, A = (e — 3p)/T 4 , is sizeable over a considerable range of 
temperatures above the deconfinement temperature, which may also indicate that the 
plasma is strongly interacting. For the ideal gas with massless particles, the interaction 
measure is zero. A nonzero value is probably due to either the interaction or the mass. 
In order to see the contribution of the particle mass to the interaction measure, in this 
paper we study a system without any interactions, i.e., an ideal gas with massive 
particles. After assembling the standard formulas of the quantum statistics with 
relativistic energy, we calculate the energy density, the pressure and the interaction 
measure. We find that their dependences on temperature reproduce the qualitative 
features of the lattice result. The interaction measure is nonzero for an ideal gas, 
which demonstrates that the particle mass contributes to the interaction measure. By 
our estimate, in the interaction measure obtained by the lattice calculation, quark 
mass contributes less than (40-50)%. There are sizeable residue interactions in the 
deconfmed phase. 
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1. Introduction 

A new state of matter, the quark gluon plasma (QGP), has been expected since the 
1970s of last century [TJ El [3]. Its experimental production in heavy ion collisions is 
later well established by the evidences at RHIC [4], [5j [T] . Many discussions have been 
focused on the question whether this plasma is weakly or strongly coupled [8]. 

The Hard Thermal Loop effective theory (HTL) is based on a physical picture of the 
QGP as a gas of weakly coupled quasi-particles — quarks or gluons with temperature 
dependent effective masses and couplings. It can well reproduce the thermodynamics 
from lattice data above 2T C , where T c is the deconfinement temperature P, [TDl EH 
[T2T [13j [H] . Even an ideal gas with thermal gluon mass can describe the gluon lattice 
data for thermodynamical quantities [15]. Thus, the thermodynamics obtained by lattice 
QCD appears to be consistent with a weak coupling picture of the QGP. However, the 
elliptic flow results at RHIC and the LHC suggest a rapid thermalization and a very 
small value of the viscosity-over-entropy-density ratio i]/s, which is inconsistent with 
the weak coupling calculations based on kinetic theory [§]. Such a small ratio can be 
obtained by the AdS/CFT correspondence at strong coupling |161 IT?] . Thus, the small 
value of i]/s supports a strongly coupled QGP (sQGP). 

It is intriguing that one system can be estimated as weakly coupled for some 
phenomena while strongly coupled for others. It is still conceivable since the coupling 
constant runs with the typical momenta exchanged in the interactions on different space- 
time scales [8]. 

The interaction measure, A = (e — 3p)/T 4 , defined in terms of the energy density 
e and the pressure p, is usually used to measure the interaction between partons in 
the lattice QCD. For an ideal gas with massless particles, energy density is equal 
to 3 times pressure and A = 0. It is thought that sizeable A may indicate strong 
interaction. The lattice QCD calculations indeed find that A is sizeable at temperature 
region T = (1 — 2)T C [18], and it appears to support a strong coupling picture for 
the QGP. However, a non-interacting quasi-particle description with only temperature- 
dependent effective mass can also well explain the behavior of the interaction measure 
from lattice [T5l [T9T 120] . The aim of the present study is to give another explanation 
on the interaction measure obtained by lattice. We study a system without any 
interactions, i.e., an ideal gas with massive particles. We first calculate the energy 
density, the pressure and the A. We find that their dependences on temperature 
reproduce the qualitative features of the lattice result. A is nonzero for an ideal gas, 
which demonstrates that the particle mass contributes to the interaction measure. We 
then further show that, in contrary to [T5l [19], ED], in our study the particle mass can 
not fully explain the interaction measure obtained by lattice. There are sizeable residue 
interactions in the deconfined phase. It may give us some new insight into the nature 
of the plasma above the deconfined point. 



The Contribution of the particle mass to the interaction measure 



3 



2. Thermodynamics of an ideal gas with massive particles 

Since both the interaction and the particle mass contribute to the interaction measure, 
in this section we will study a system without any interactions, i.e., an ideal gas with 
massive particles, to separate their contributions. The reason we use an ideal gas is 
that we want to see the pure contribution of the particle mass when there are not any 
interactions. 

For an ideal gas with particle mass m at an equilibrium state described by (T, /i, V), 
the standard formulas of the quantum statistics give the energy density 

/A 3 
(>}3^ F ' B ' ^ 

and the pressure 

where g is the degeneracy factor, and /f,b stands for the single-particle distribution 
function with 

f F (s,/i,T) = e(£ _ M)/T + 1 for fermions, (3) 
f B (e,ii,T) = „ (£ _ M)/r _ 1 for bosons. (4) 



c 



Using the relativistic dispersion relation e = \J\~fz \ 2 + m 2 , we obtain 

d 3 k e 2 — m 2 „ r d 3 k , m 2 



3P = 9 J (2vr)3 S /f ' b = 9 J (2^ (£ " V )/f ' B ~ " (5) 
The above equation tells us that the presence of mass always reduces the pressure 
below e/3. It is easy to understand because massive particles move slowly at a given 
temperature and thus the thermal pressure reduces. The interaction measure is then 

/d^ Iv m 2 

In the following we fix the chemical potential /i — for comparison with lattice 
results. At n — 0, the pressure in equation (2), the energy density in (1) and the 
interaction measure in (6) are only functions of T and m. In this paper the integrals for 
the quantum statistics are done numerically. 

Besides, we also do a Boltzmann statistics, i.e., the quantum distributions are 
replaced by the Boltzmann distribution in the above equations, which can be analytically 
done. For the pressure, a technique of integrating by parts with the relation J^e~ £//T = 
— i|e~ £//T is used. Finally we obtain 

c\n 3T 4 T 4 

e(T;m) = T^- P (T) = ^K 2 {x) + ^^(x), (8) 

e-3p=^x 3 K 1 (x), (9) 
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Figure 1. Energy density, pressure and interaction measure are plotted, (a) the 
lattice QCD result of strongly interacting matter from figure 7 in reference [TS] with 
m u( i « 0.1m s . Diamonds, squares and circles are for energy density, pressure and 
interaction measure, respectively. (Reprinted figure with permission from Bernard C 
et al, Phys. Rev. D, 75, 094505, 2007. Copyright (2007) by the American Physical 
Society.) (b) the result of an ideal gas with particle mass m = l.OGeV. Red lines 
for bosons, blue lines for fermions and green lines for Boltzmann approximation. The 
arrows near the right axis indicate the corresponding Stefan-Boltzmann limits. (Color 
online) 

with x = m/T. 

At fixed mass, e.g. m = 1.0 GeV, p{T), e(T) and A(T) are shown in figure 
1(b), for both quantum statistics and Boltzmann statistics. For easy of comparing, 
the lattice results are shown in figure 1(a). The lattice studies the equation of state 
of a real QCD matter at // = 0. In figure 1(a), a continuous increase of the energy- 
density and the pressure is seen, indicating a crossover between the hadron gas and 
the quark gluon plasma, which is confirmed by later reference [2TJ. There are more 
recent results [22 | [23 | 124"] on the QCD equations of state than shown in figure 1(a) from 
reference [18] . In the temperature region T c < T < 2T C (T c = 170MeV), the interaction 
measure, denoted as / in figure 1(a), is sizeable and does not vanish. 

The analysis is as follow: 

(i) In figure 1(b), e/T 4 and p/T A are increasing with T and saturate at high 
temperature. At high temperature limit, p and e approach Stefan-Boltzmann 
limits marked as arrows in the figure. (The SB limits are given in Appendix A.) 
A/T 4 is nonzero and peaks at the temperature region where the pressure and 
the energy density increase rapidly, which is similar to the lattice result. It is 
interesting that, even an ideal gas has the similar trend as the lattice results for the 
QGP thermodynamics. Besides, another kind of a more complicated weak-coupling 
calculations from the kinetic theory[9j QUI EU [121 E2 EH] can describe the lattice 
data not only the trend but also quantitatively. In addition, an ideal gas with 
temperature-dependent mass also leads to good agreement with lattice [I5j HHJ EH] . 
The system used here is quite simple, which only reproduce the qualitative features 
of thermodynamic quantities. The aim of this paper is to investigate how much 
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Figure 2. The mass dependence of the interaction measure. Red, green, blue lines 
are for mass 1.0, 0.8, 0.2 GeV respectively. This plot is calculated from Fermi-Dirac 
statistics. (Color online) 

the particle mass contributes to the interaction measure. By comparing our results 
with the lattice results, we can obtain the pure contribution of the particle mass. 
A quantitative analysis will be given in the subsequent paragraphs. 

(ii) As discussed before, for an ideal gas with massless particles, A = 0. In figure 1(b), 
it can be seen that A is nonzero for an ideal gas with massive particles. The only 
difference between the two systems is the particle mass. Thus, the lesson we learn 
from here is that, a nonzero A does not certainly mean strong coupling since the 
particle mass contributes to this quantity. 

(iii) Here we use m — 1.0 GeV as an example to show the behavior of the interaction 
measure. Similar to the lattice result, there is a peak at the temperature region 
where the energy density and the pressure grow fast. By changing the particle 
mass, the curve of the interaction measure will shift its position, as shown in figure 
2. The A peak shifts to low T for small mass and to high T for large mass. And 
the peak height is independent on mass for an ideal gas. This statement is valid 
for an ideal gas with either quantum statistics or Boltzmann statistics. 

(iv) In figure 1(b), the peak height of A/T 4 for the ideal gas is 0.057 (this value is 
for Fermi-Dirac statistics, 0.059 for Boltzmann statistics, 0.061 for Bose-Einstein 
statistics), which is much smaller than the lattice calculation. In figure 1(a), the 
peak height is about 6.5, while newer results from HotQCD collaboration favor a 
value of about 5.0 for the peak height [22], and the continuum extrapolated result 
of the Wuppertal-Budapest collaboration indicates a peak height of about 4.1 [23] . 
The great disparity in peak height between figure 1(a) and figure 1(b) can be 
partly attributed to the number of degrees of freedom of constituent particles. The 
degeneracy factor g = 1 in the ideal gas system in figure 1(b). However, g in the 
QCD matter considered by lattice refers to a large number of degrees of freedom, 
including gluon and three flavor quark, anti-quark which further count color, spin 
degeneracy. 

Taking the pressure for an example. In an ideal multicomponent gas, the partial 
pressure of one component is d x p m , where the particle degeneracy factor d 
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represents the number of the internal degrees of freedom for the particle with mass 
m, p m represents the pressure without any internal degrees of freedom. The total 
pressure of an ideal multicomponent gas can be expressed as 

p(T; mi, ... , m n ) = d t p(T, mi) + . . . + d n p{T, m n ). (10) 

Similarly, 

e(T;mi, . . . ,m n ) = die(T,mi) + . . . + d n e(T,m n ), (11) 

and the interaction measure is 

A(T; mi, ... , m n ) = d t A(T, mi) + . . . + d n A(T, m n ). (12) 

In lattice calculations, only quark mass contributes to the interaction measure since 
gluon is massless, according to equation (6). For three flavor QCD, there are two 
masses: m u d and m s . A for these two masses are peaks at different positions. 
According to (12), A for the mixed gas will be a double-peak curve. 
The maximum value of A in equation (12), denoting as A max , is less than 

A max < d^ZT + • • • + d n^ZT, (13) 

where A™ ax represents the peak height of the A curve with fixed mass m. Since the 
peak height is independent on mass for an ideal gas, then A™ ax = A™ ax = Ao = 

0. 057. Thus we can infer that 

A max < d t otalA Q , (14) 

where d to tai — d q + d q = 36 for three flavor QCD, with the quark degeneracy 
d q = dq = N c N s Nf = 18 [25]. Finally we get 

A max < rf total A = 2.1. (15) 

The value at the right hand side of the above equation is about 40% of the lattice 
result 5.0 reported by HotQCD Collaboration [22], and about 50% of the lattice 
result 4.1 reported by Wuppertal-Budapest Collaboration [23]. Therefore, from our 
analysis we can conclude that the pure contribution of the particle mass is not 
larger than (40-50)%. 

In this paper, we only focus on the peak height which depends on the interaction. 
The peak position may also shift by additional interaction. Thus the peak position, 

1. e., the temperature where the interaction measure gets maximum, for an ideal gas 
can not be used to compare with that of the lattice QCD. 

3. Conclusions 

We use an ideal gas with massive particles to calculate the energy density, the pressure 
and the interaction measure, and compare them with that of the lattice results. An 
ideal gas has the similar trend as the lattice results for the QGP thermodynamics. 
We reproduce the qualitative features of the lattice results on the energy density, the 
pressure and the interaction measure by an ideal gas with massive particles. The 
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interaction measure is nonzero for an ideal gas, which demonstrates that a nonzero 
A does not certainly mean strong coupling. Particle mass can make the interaction 
measure nonzero. The A measured in the lattice QCD includes the contribution from 
the particle mass because they indeed use quark mass in the lattice calculation [18j [26] . 
After counting the degeneracy number of the QGP, A contributed by the particle mass 
explains less than (40-50)% of the lattice result. We infer that the other contribution 
comes from interaction. That means nonzero A from lattice calculation does not mean 
strong coupling but does not exclude a strong coupling picture of QGP. It is still hard to 
know how strong is the coupling since it is hard to know the relation of the interaction 
measure with the coupling. The question how strong is the coupling needs further study. 
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Appendix A. Stefan-Boltzmann limits 

Stefan-Boltzmann limit is the value for a non-interacting gas with massless particles. 
For Boltzmann particles, since x 2 K 2 (x) — )■ 2 for x — > 0, equation (7) (8) becomes 

p(T; m) = ?—x 2 K 2 (x) -> ^, (A.l) 



27T 2 V ' 7T 2 

,dp ^ , 3T 4 

' V 2 



e(T-m)=T^-p(T)^—. (A.2) 



Therefore, 



— > — « 0.303 for Boltzmann particles, (A. 3) 

T 4 7T 2 



as marked in figure 1(b). 

For quantum particles, ref [25] gives 



e 7tt 2 

—j — > « 0.288 for fermions, (A. 4) 

T 4 8 30 K J 

e tt 2 

—r ->■ — w 0.329 for bosons, (A. 5) 

T 4 30 y ' 



which are also marked in figure 1(b). 
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